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We show that the synchronization transition of a large number of noisy coupled oscillators is an 
example for a dynamic critical point far from thermodynamic equilibrium. The universal behaviors 
of such critical oscillators, arranged on a lattice in a d-dimensional space and coupled by nearest 
neighbors interactions, can be studied using field theoretical methods. The field theory associated 
with the critical point of a homogeneous oscillatory instability (or Hopf bifurcation of coupled oscil- 
lators) is the complex Ginzburg-Landau equation with additive noise. We perform a perturbative 
renormalization group (RG) study in a 4 — e dimensional space. We develop an RG scheme that elim- 
inates the phase and frequency of the oscillations using a scale-dependent oscillating reference frame. 
Within a Callan-Symanzik RG scheme to two-loop order in perturbation theory, we find that the 
RG fixed point is formally related to the one of the model A dynamics of the real Ginzburg-Landau 
theory with an 0(2) symmetry of the order parameter. Therefore, the dominant critical exponents 
for coupled oscillators are the same as for this equilibrium field theory. This formal connection with 
an equilibrium critical point imposes a relation between the correlation and response functions of 
coupled oscillators in the critical regime. Since the system operates far from thermodynamic equi- 
librium, a strong violation of the fluctuation-dissipation relation occurs and is characterized by a 
universal divergence of an effective temperature. The formal relation between critical oscillators and 
equilibrium critical points suggests that long-range phase order exists in critical oscillators above 
two dimensions. 

I. INTRODUCTION 

Equilibrium systems consisting of a large number of degrees of freedom exhibit phase transitions as a consequence 
of the collective behavior of many components [USE]- The universal behaviors near critical points have been studied 
extensively using field theoretical methods and renormalization group (RG) techniques [3| r5|. f6|. fTt f8| fQ^ fTOl fTTl IT21 IT3] . 
In systems driven far from thermodynamic equilibrium, collective behaviors can lead to dynamic instabilities and non- 
equilibrium phase transitions |14l 1151 [TBI 117) . While the study of non-equilibrium critical points has remained a big 
challenge, RG methods have in some cases been applied [T51 ITUl I2TJ1 |2"T1 |2"2"] . 

An important example for non-equilibrium critical behavior is a homogeneous oscillatory instability or Hopf bifur- 
cation of coupled oscillators [23 . Such instabilities are important in many physical, chemical, and biological systems 
[2~4"ll25| . From the point of view of statistical physics, phase coherent oscillations result as the collective behavior of a 
large number of degrees of freedom in the thermodynamic limit. For a system of finite size, fluctuations destroy the 
phase coherence of the oscillations and the singular behaviors characteristic of a Hopf bifurcation are concealed. 

This general idea can be illustrated by individual oscillators arranged on a lattice in a (i-dimensional space and 
coupled to their nearest neighbors. As a consequence of fluctuations in the system, each oscillator is subject to a 
noise source. For small coupling strength and as a result of fluctuations, the oscillators have individual phases and 
exhibit a limited coherence time of oscillations. In the thermodynamic limit, a global phase emerges beyond a critical 
coupling strength where oscillations become coherent over large distances. As this critical point is approached from 
the disordered phase, a correlation length, which corresponds to the characteristic size of the domains of synchronized 
oscillators, diverges. At the same time, oscillations become coherent over long periods of time and long range phase 
order appears. The divergence of the correlation length allows a description of the critical behaviors of spatially 
extended systems in terms of continuous field theories and is characteristic of scale-invariance in the critical regime. 
It constitutes the foundation of the renormalization group (RG) theory, which explains the emergence of universality 
at critical points H EJ [TTJ US US] . 

In this paper, we perform an RG study of the critical behaviors of a collection of oscillators, distributed on a d- 
dimensional lattice and coupled via nearest neighbors interactions. We approach the critical point from the disordered 
phase. Close to criticality, the large scale properties of the array of oscillators are described by a dynamic field theory 
that is given by the complex Ginzburg-Landau equation with an additive noise term. We apply field theoretical 
perturbation theory in a d = 4 — e dimensional space and introduce an RG scheme that is appropriate for the study 
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of coupled oscillators. 

The outline of the paper is as follows: In Section [TTJ we present the general field theoretical framework for the 
complex Ginzburg-Landau equation. We introduce in Section |III| an oscillating reference frame that is essential to 
define the RG procedure for oscillating critical systems and in which the phase and frequency of the oscillations are 
eliminated. The correlation and response functions of critical oscillators in mean field theory arc discussed in Section 
|IV| These mean field results are relevant above the critical dimension d c — 4. For d < 4, mean field theory breaks down. 
We discuss in Section [V] the renormalization group of the complex Ginzburg-Landau field theory using a Wilson's RG 
scheme for which the renormalization procedure for critical oscillators can be introduced most clearly. One-loop order 
calculations in perturbation theory are presented, but further calculations are necessary to characterize the correct 
qualitative structure of the RG flow. In Section |VI[ we present the Callan-Symanzik's RG scheme for systems of 
coupled oscillators and calculate its beta functions as well as its complete RG flow and fixed points to two-loop order 
in perturbation theory. The physical properties of the oscillating system are characterized by correlation and response 
functions. In Section fVlI| we discuss the asymptotic behaviors of these functions in the critical regime using the RG 
flow and applying a matching procedure. The formal relation of the RG fixed point for critical oscillators to the fixed 
point of a real Ginzburg-Landau theory (which satisfies a fluctuation-dissipation (FD) relation) leads to an emergent 
symmetry at the critical point. However, the system operates far from thermodynamic equilibrium and breaks the 
FD relation. The degree of this violation can be characterized by the introduction of a frequency-dependent effective 
temperature which diverges with a universal anomalous power-law at the critical point. We conclude our presentation 
with a discussion of the general properties of critical oscillators, their relations to equilibrium critical points and 
possible experimental systems for which the critical behaviors discussed here could be observed in the future. 

II. FIELD THEORY OF COUPLED OSCILLATORS 

A. Complex Ginzburg-Landau field theory 

The generic behavior of a nonlinear oscillator in the vicinity of a Hopf bifurcation can be described by a dynamic 
equation for a complex variable Z characterizing the phase and amplitude of the oscillations [24 . This variable can 
be chosen such that its real part is, to linear order, related to a physical observable, e.g. the displacement X(t) 
generated by a mechanical oscillator: X(t) — Re(Z(i))+nonlinear terms. In the presence of a periodic stimulus force 
F(t) = Fe~ lult with a frequency u> close to the oscillation frequency ojq at the bifurcation, the generic dynamics obeys 

d t Z = -(r + iuj Q )Z - (u + iu a ) \Z\ 2 Z + A^e^Fit). (1) 

For F — and r > 0, the static state Z = is stable. The system undergoes a Hopf bifurcation at r — and exhibits 
spontaneous oscillations for r < 0. The nonlinear term, characterized by the coefficients u and u ai stabilizes the 
oscillation amplitude for u > 0. The external stimulus appears linearly in this equation and couples in general with 
a phase shift 9 [3Hl|57j. In the case of a mechanical oscillator, the coefficient A has units of a friction. 

Coupling many oscillators in a field theoretic continuum limit leads to the complex Ginzburg-Landau equation 
[2"5l |2"5] with additive noise and external forcing terms: 

d t Z = -(r + iuj )Z + (c + ic a )AZ - (u + iu a ) \Z\ 2 Z + h^e^F + 77. (2) 

Here, the complex variable Z(x, t) becomes a field defined at positions x in a <i-dimensional space and A denotes the 
Laplace operator in this space. The coefficients c and c a characterize the local coupling of oscillators and the effects of 
fluctuations are described via a complex random forcing term ?y(x,t), which will be chosen Gaussian with zero mean 
value, i.e. (?7(x, t)) = 0. As far as long time and long wavelength properties are concerned, the correlation times of 
the noise can be neglected and white noise can be used. 

For a vanishing external field F(x,t) and in the absence of fluctuations, Eq. ^ is invariant with respect to phase 
changes of the oscillations: 

Z Ze^. (3) 

This symmetry reflects the fact that only phase-invariant terms contribute to the dominant critical behaviors studied 
here. Indeed, the Hopf bifurcation is associated with the emergence of a non-zero oscillatory mode, which dominates 
the critical behaviors and for which time-translational invariance and phase invariance are equivalent. The noise 
correlations can therefore be chosen such that they respect phase invariance in the problem: 



(77(x,t)^(x',i')> = 

(7?(x,t)?7*(x',t')) = 4ZW d (x-x')<S(i-i'). 



(4) 
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Here D is a real and positive coefficient characterizing the amplitude of the noise, and 5 and 8 d represent Dirac 
distributions respectively in f and d dimensions. 



B. Physical correlation and response functions 



Since the physical variables of interest are real, we decompose the complex fields Z and F into their real and 
imaginary parts by Z — tpi + iip2 and F = F\ + iFi. We focus on two important functions which characterize the 
behavior of the system, namely the two-point auto-correlation function C a p and the linear response function Xa/3 to 
an applied external forcing term. They are defined as 

C Q/3 (x - x', t - t') = (^ a (x, i)^(x', t')) c 

(Vfete *)) = / d d x' dt' Xq/3 (x - x', t - t')Fp(x.', t') + 0(\F\ 2 ), (5) 

where (...) c denotes a connected correlation function j51j. Because of phase invariance, these functions obey the 
following symmetry relations: 

Cn = C 22 and C 2 i = -C 12 , (6) 
with similar relations for the linear response function Xafi- 



C. Field theoretical representation 



The correlation and response functions C a p and \af3 can be conveniently expressed using a field-theoretical for- 
malism. We introduce the Martin-Siggia-Rose response field ip a j3U] and apply the Janssen-De Dominicis formalism 
[3~Tl [35] to write the following generating functional: 



Z f 



lot i 1 a 



V [i/, a ] V 



exp < Sp 



d d xdt 



where the "physical" action Sp is given by[ 



1pa,i>a 



d d xdt ^Dlp a ljj a - 1p a [d t tp a + (R a p +WoEap)4>P + Uapl/jplp^lpy]} 



(7) 



(8) 



Here R a p = (r — cA)S a/ 3 — c a A e a p and U a p = uS a p + u a e a p, with £21 = —£12 = 1 and e\i — £22 = 0. The field I a is 
related to the external force F a in the equation of motion (Eq. by 

(9) 



where £l a p(9) denotes the rotation matrix by an angle 9 in two dimensions: 



Qa/3\9) — 



cos(6>) -sin(0) 



sin(0) cos(0) 

Correlation and response functions are given by derivatives of the generating functional. We have 

C Q/3 (x-x',t-i') 



SI a (x,t)6I (x>,t 



1 H\ II a ,I<*=0 



Xa p(x-x',t-t') = A-^pie) 



S 2 \nZ 



<5/ Q (x,i)<5/ 7 (x',i') 



\i„,i a =o- 



(10) 

(11) 
(12) 



III. FIELD THEORY IN AN OSCILLATING REFERENCE FRAME 



A. Amplitude equation 



The frequency ui Q and the phase 9 can be eliminated from Eq. ([2]) by a time-dependent variable transformation of 
the form Y = e luJot Z, H = e luJ ° t A~ 1 e F and C = e lu>at r], where Y denotes the oscillation amplitude, H is a forcing 
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amplitude and ( is a transformed noise that has the same correlators as rj. This procedure leads to the amplitude 
equation 



d t Y = -tY + (c + ic a ) AY - (u + iu a ) \Y\ 2 Y + H + (. 



Defining two real fields 4> a by Y — <f>i + i<p2, Eq. (13 1 reads 



9t4>a = -Rap4>p - U a p(j)p4>~<4>~1 + H a + £ a , 



(13) 



(14) 



where H = H x + iH 2 . 

The correlation and response functions G a p — (4> a 4>p} c and 7 Q /3 = (<Pa4'p}c of the fields 4> a are related to the 
physical correlation and response functions by 



C a/3 (x,i) = tt aa (~uj t) G CT/3 (x,i) 

Xa/3(x,t) = A -1 tt aa (6 - UJ t)-f al3 (x,t). 



(15) 



Similar time-dependent transformations exist between higher order correlation and response functions of the physical 
fields and those calculated in the oscillating reference frame. 



B. Analogy with an equilibrium critical point in a particular case 



For the particular case where c a = and u a = 0, Eq. (14 1 becomes identical to the model A dynamics of a real 



Ginzburg-Landau field theory with an 0(2) symmetry of the order parameter [3]. The critical behavior of this theory 
at thermodynamic equilibrium has been extensively studied |33l 1341 135) . This leads, in this particular case, to a 
formal analogy between an equilibrium phase transition and a Hopf bifurcation. The amplitude Y plays the role 
of the order parameter of the transition. The disordered phase with (Y) = corresponds to noisy oscillators that 
are not in synchrony, while nonzero order (Y) implies the existence of a global phase and amplitude of synchronous 
oscillations with frequency luq. The correlation lengths and times of the equilibrium field theory correspond to lengths 
and times over which oscillators are in synchrony. The correlation and response functions C Qi g and Xap for this case 
can be obtained from those of the equilibrium field theory by using Eq. (15). Since the 0(2) symmetric theory at 



thermodynamic equilibrium obeys an FD relation, a generic relation between the correlation and response functions 
C a p and Xa/3 appears. 



C. Generating functional 

The functions G a p and r y a p 1 as well as higher order correlation and response functions, can be formally calculated 
using field theoretical techniques (see e.g. [2j [13]) . The generating functional of the theory is given by 



J at J a 



v [<t> a ] V 



exp < S 



d d xdt 



where we have introduced the Martin-Siggia-Rose response field <p a [3D]. The associated action reads: 



= d 



s 

Correlation and response functions are given by 

G a/3 (x-x',i-i') 
j a p(x.-x',t-t') 

The effective action of the theory is defined by 



l d X(lt $yD4> a <}) a - 4> a [d t 4>a + Rap4>p] ~ ^a^^a^^^} 



5 2 In Z 



5 2 \uZ 

5J a (x 7 t)5Jp(x> 7 t>) \ J -J°=°- 



r $Q,,<i>Q, = / d d xdt J a $ a + J a $ a -InZ J a ,J a 



(16) 
(17) 

(18) 
(19) 

(20) 
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with 



i / \ ShiZ ~ , . ShiZ 

$ a (x,f) = and $ a (x,<) = 

dJ Q (x,i) 6J a (x,t) 



(21) 



In order to perform calculations perturbatively, we split the action into a harmonic or "Gaussian" part, and a 
quartic or "interaction" part as: S = So + S- m t. So is given by 



S = I {} a (k) (D4> a (-k) - [iujS afj + fl Q/3 (-fc)]<^(-fc)) } 

1 r Lt 



2 j/ a m a ,(-^k), 



where 



(22) 



(23) 



I ,_ /.,/ -2D5 a/ 3 iuj8 a/3 + R al3 (-k) 







(24) 



and R a p(k) = (r + cq 2 )5 a p + c a f^s a p- In these expressions and in the following, we label k = (q,w), J fc = J ( 
denotes J r^ya > an d we use the following convention for Fourier-transforms: 



q,u; 



/(x,t) = / /(q,c)e^ qx -^ = / f(k)e lk - x . 

J q.cj J k 



(25) 



The interaction term of the action takes the form 



Graphic representations of the basic diagrams of the perturbation theory are given in Appendix |Aj 



(26) 



IV. MEAN FIELD THEORY 



Dimensional analysis reveals that for d > 4 mean field theory applies. In this case, the mean field approximation 
allows us to calculate valid asymptotic expressions for the effective action of the theory, as well as the two-point 
correlation and response functions. In the framework of the Janssen-De Dominicis formalism for dynamic field 



theoretical models, this approximation consists in substituting the saddle-node value of the path-integral (16 1 to the 
full functional generator Z. We obtain 



7inf 



J a ,J a =exp^S 



d d xdt [J a tf? + J a 



where S is given by ( 17 1 and 0™ and </>™ satisfy the stationarity conditions at the saddle node 



SS 



- J a and 



SS 



(27) 



(28) 



which give the mean field dynamic equations 

Ja = -d t 4>a f + Rpaftp* + U paftjj* (j>™ <f>™ 1 + 21]^^^^ 



(29) 



The correlation and response functions can be obtained most easily by first eliminating the nonphysical field cf>™ f and 
writing a mean field equation for 0™ f only: 



2DJ a = [-d t 6 a/s + Rpc + Up^flf + 2U 0a ^f^ 



(30) 
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The mean field generating functional for the field 0™ f obeying Eq. ( 30 1 is obtained from ( 27 1 by eliminating the field 
6™*. It reads: 



J a j J a 



' x:p { ~ Id J ddxdt ( 9i ^" lf - F <* -Jo) 2 + 4C f } 



where 



Note that for J a = the stationary condition of the generating functional (31 1 leads to 

5 t C f = ~R a ^f - Uaf,<l>tf41?fif + J a , 



(31) 

(32) 
(33) 



which is compatible with Eq. ( 30 ) 



The linear response and correlation functions are obtained as 



G^(x-x';t-i') = 



8j {x>,t>y J °< J "= ' 



(34) 



Calculating these functions and applying the time-dependent transformations (15), we obtain the physical correlation 
and response functions in mean field theory 



f 1 

X^g(q,w) = j^{{(-iuj + R) cos 9 + fl sin 8] 5 a/ 3 + {(-iid + R) sin 9 - fl cos 6) e a/ 3} 



2D (u? + ttl + R 2 2iton 

-2iu,tt u} 2 + fl 2 ) + R 2 ' 



(35) 



where ii = r + cq 2 , SIq = + c a<l 2 and A = (— iu> + R) 2 + Qg. The diagonal elements of these matrices are given by 



xn(q^) 

Ci m i f (q^) 
and the non-diagonal elements by 

cr 2 f (q,c) 



1 

2A 



R — i(u — Cl ) R-i(u + £l ) 
D D 



R 2 + (Lu-n Q ) 2 R 2 + (to + n y 



(36) 



i 

2A 



R-i(u-fl ) R-i(uj + Q ) 
iD iD 



R 2 + (cj-n Q ) 2 R 2 + (aj + n ) 2 ' 



(37) 



Note finally that for r — and for the critical mode (q = 0, w = 0), the response function of the system is nonlinear 
even at small amplitudes. We have indeed in this case 



S(\Y(O,0)\) ml oc \H\ 1/3 . 



(38) 



V. 



WILSON'S RENORMALIZATION SCHEME 



For d < 4, mean field theory breaks down and another approach is necessary to investigate the critical behaviors of 
the theory. We apply perturbative renormalization group (RG) methods using an e expansion near the upper critical 
dimension d c = 4 (d = 4 — e) [7, 36, 37]. We present here the RG structure of the theory within a Wilson's momentum 
shell RG scheme adapted to the renormalization of the complex Ginzburg-Landau field theory. This scheme has the 
advantage to be conceptually transparent and to provide a clear physical interpretation to the calculations. However, 
for calculations beyond one-loop order, this technique is less suited than the Callan-Symanzik RG scheme. The 



adaptation of the latter to the complex Ginzburg-Landau theory is presented in Section VI 
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A. Renormalized fields 



The renormalization procedure within the Wilson's scheme is performed as follows. We start from a dynamic 
functional of the theory with a small distance cut-off A in the integrals over wave vectors. This cutoff corresponds to 
an underlying lattice of mesh size a ~ 2w/A. We interpret this as a microscopic theory with an action of the form 
^ , and associated quantities are labeled with an superscript "0" . We calculate the effective action in an oscillating 
reference frame at a given scale A/6, where b — e l is a dilatation coefficient larger than 1. This reference frame is 
defined such that, described in terms of effective or renormalized parameters, the effective action has the structure 



(17 1. The renormalized quantities that satisfy this requirement can be expressed as 



q = bq° 

uj = tfZ^fyuj 

o (q,t) = b-™ y%(6)z w (&)tt Q/3 (^ (^v>°(q , t°) 

Mq,t) = ^Z^b)Z u (b)Q aP (9{b) +w (&)t)^(qV )- (39) 

Here we have introduced scale-dependent Z factors for the renormalization of the two dynamic fields (2^(6) and 
Z^{b)) and for the frequencies {Z^ib)). In addition to the usual RG transformations and scale dilatation q = 
bq(°\ the complex Ginzburg-Landau theory requires us to perform time-dependent transformations between bare 
and renormalized fields. Indeed, the effective theory is described in a reference frame that oscillates with effective 
frequency Ojo(b) and phase 6(b) relatively to the reference frame of the bare theory. The definition of the fields cf> a 



and (j) a takes this relative rotation into account by terms involving the rotation matrix f2 Qj g (Eq. (lOl). The scale- 
dependent oscillating reference frame represents a key element of the RG procedure for oscillating systems. For the 
microscopic theory, cD (l) = oj®, 9(1) = and the fields ip^ and -0° coincide with the fields introduced in Section II. 

B. Renormalization group flow 

In order to determine the behavior of the effective parameters under renormalization, we determine their variations 
with respect to a small changes of the dilatation Sb/b (or 81). Integrating over the momentum shell of wave vectors 
in the interval [A/e , A], we obtain an effective action that, before rescaling, reads: 



U Q >J = / / {(D + 8D)^ a -^ a [i(l + 8X 1 )oj+(r + 8r) + (c+8c)q 2 ]^ a 

-(j> a [1SX2UJ + Sujq + (c a + Sc a )q 2 ] e a p<j>p \ - / (U a(i + 8U a p)<f> a <l>p^i<t>i- 

> -'kikaks 



(40) 



Here, we have introduced variations of the parameters under this procedure, which can be calculated perturbatively 
from the "one-particle irreducible" (1-PI) diagrams of the theory. Other terms are either forbidden by symmetry 
properties of the theory or irrelevant in the infra-red (IR) limit. Two new terms appear, which were absent from Eq. 



(17 1 and which correspond to the coefficients SX2 and 5ujq. They reflect the renormalization of frequency and phase. 



These terms are absorbed by time-dependent variable transformations of the fields: 

<#*(q>*) = n «/8 (8&bt) 4>f}(q,t) 

&(q,t) - na/j(^ * + ^)^/s(q.*). ( 41 ) 

where 

66 = arctan [-<5A 2 (1 + (SAi)" 1 ] 
Su)^ = Sluq cos S9 + (r + Sr) sin S9. (42) 

Rewriting the effective action in terms of the new fields <f) l a and <p l a requires a redefinition of the changes of all 
parameters. The so defined parameter changes are label by a superscript "i" in the following. We finally rescale 
all lengths by q' = e 5l q, and introduce three Z factors at the scale I + SI, such that the effective action retains its 
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form (17) under an RG step. Furthermore, we impose that D(l + SI) — D(l) and c(l + 51) = c(l), i.e. that the 
coefficients c and D remain constants under renormalization. As a result, the variations of the parameters under a 
small renormalization step are given by 



c a (l + 5l 
r(l + 6l 
U a/3 (l + 5l 
Z u (l + 5l 
Z^l + Sl 
ZM + Sl 



= [c a (l) + 6d a ] (1 + Sd/c)- 1 
= e 251 [r(0+<5r>] (l + Sd/c)' 1 

= e^ sl [U afi (l) + SUlp] (1 + Sd/c)- 2 ^ + 5D l /D)(l + SX i y 

= Z^il + Sj/cy^l + SX) 

= Z(l)(l + 5d/cf(l + SD'/D)- 1 ^ + SX)- 1 

= Z(l)(l + 8&/c)(l + 5D { /D)(l + 5X)-\ 



(43) 



where 



5D 1 = 5D 

5X l = (G - 1) + JAiC - 5X 2 S 

5r l = r(C-l)+5rC-(w + 6u )S 

5d = c(C-l)+5cC-(c a + 5c a )S 

5c) a = c a (C-l) + 5c a C+(c+5c)S 

sulp = u aa (n aP (56)-i aP ) + su a M*p(se). 



(44) 



Here, C and 5* denote cos 56 and sin (56*, respectively. The evolution of the effective frequency and phase that enter 
the transformations ( 39 1 are given by 



6(1 + SI) = 6(1) + 56 
u (l + 5l) = e 2Sl [w (Z) + <foj] (l + Sj/cy^l + SX). 



(45) 



C. Correlation and response functions 

Using the RG transformations of the effective parameters, as well the correlation and response functions of the 
renormalized fields <f> a and <j) a , we can write expressions for the physical correlation and response functions, where 
the transformations (39 1 have been used: 

C a0 (q,t) = 6 2 Z (6)- 1 Z w (&)- 2 r! QC7 (-^ o (fo)i)G« /3 (6q, 6- 2 Z w (6)-4) 

Xa/s(q,*) = y/z4b)Z^b)~' V Z^(bY 2 fl aa (6(b) - tu (b) t) 7 * (6q, b~ 2 Z a {b)~H). (46) 

Here, the superscript "i?" indicates that the functions have to be calculated using the renormalized set of parameters. 
The scale-dependent renormalized frequency is given by 

wd(6) = b- 2 Z LJ (b)- 1 oj (b), (47) 

and 6(b) describes the scale-dependent phase lag between external forcing and response of the system. 

At the fixed point of the RG, the theory is scale invariant and, consequently, the Z factors exhibit simple scaling 
relations as a function of b: 

Z u (b) = b z - 2 
Z^b) = b- 2 ^ +r ' 

Z^(b) = &-2(*-2)+« . (48) 

These relations define three independent critical exponents of the theory, z, r\ and fj. A further critical exponent v is 
associated with the positive eigenvalue of the linearized RG equations around the fixed point. 
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D. Renormalization of the time and independent critical exponents 

In dynamic RG procedures, there is in general a freedom to choose some parameters constant while others are 
renormalized in a non-trivial way. In the system described here, we choose to renormalize the time (i.e. the frequency 
coordinate uj in Fourier space) and to keep the parameters c and D invariant under renormalization. In other 
cases, different choices are commonly used. For example, renormalizing the dynamic model A with 0(2)-symmetry 
at thermodynamic equilibrium, the coefficient D is usually chosen to change under renormalization, while time is 
simply rescaled [531 133 . This model corresponds to the particular case of the amplitude equation of the complex 
Ginzburg-Landau theory where both c a and u a are equal to zero. It is described by a dynamics of the form 

d t <i> a = -D^ + C a , (49) 

and relaxes towards a thermodynamic equilibrium. Note, that since this model satisfies an FD relation, the noise 
strength D appears as the mobility coefficient in the dynamics. Eq. f49| shows that both choices, renormalizing the 
time or renormalizing D 7 are equivalent. The complex Ginzburg-Landau theory discussed here does not obey an FD 



relation. Therefore, a factorization of the coefficient D as it is done in Eq. (49 1 would be artificial in this case, and 
would enforce to redefine all other parameters. Without this factorization, the structure of the theory imposes to 
renormalize the time if the coefficient c is kept constant. 

Note also that the absence of an FD relation in the theory changes the structure of the RG equations as compared 
to an equilibrium 0{2) model. Indeed, the FD relation imposes a constraint on the renormalized quantities of the 
equilibrium theory, which can be written as 

Zuy/Z^Zf- EE 1, (50) 

and which implies the following relation between the critical exponents [13 : 

z=2+(fj-r])/2. (51) 

In the non-equilibrium case considered here, such a constraint does not exist and four truly independent critical 
exponents are present in the theory as compared with three in the 0(2) dynamic model A. 



E. Results to one-loop order 



To one-loop order in perturbation theory, the Z factors, as well as the phase factor and the parameter c a , are not 
renormalized. We define the following reduced parameters 

f = r/cA 2 ; w = w /cA 2 ; c a = c a /c 



u = u/(cA 2 ) 2 ; u a = u a /{cK 2 ) 2 ; 
ADA d 

V 



(47r) d / 2 r(d/2) ' 

to express the RG flow in d = 4 — e. It is given by three coupled equations 



dr 
~dJ 
du 

~di 

du a 
~dJ 



= 2f+2V- 



= eu-V 



1 + f 

2 „-,2 



(u 2 - uj)(l + r) + 2uu a c a 

(1 + f) K + (l+f)2] 



m 2 



(i + fy 



v- 



i 



(u 2 - u 2 a )(\ + r) + 2uu a Cg 
(1 + f) [cg + (l+ f )2] 



- &D 



(1 + f) 2 



and a fourth one associated with the renormalization of the oscillation frequency: 

dujQ 



dl 



2V- 



1 + f 



(52) 



(53) 



(54) 



This last equation has to be integrated after the system (53 1 has been solved. 
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Since the coefficient c a is not renormalized to first order, we find one infrared stable fixed point for each of its 
values, ft is given by: 



e e a -* - e 

■ — i u — — - and u„ — c„ — — . 

5' 52? 52? 



(55) 



Writing r — r* + <5r, u — u* + Su and u Q = u* + 5u a , the linearized RG equations at the fixed point are given by: 



d 
Jl 



Sr \ /2(l-e/5) 22? 
5u = -e 

Su a \ -4c a e/5 -el 







) 


Su 







(56) 



The RG flow of the theory is three dimensional. We show in Fig. [T|the qualitative RG flow projected on the plane 
(r,u). Finally the critical exponents to the one-loop order read: 




FIG. 1: Qualitative representation of RG flow to one-loop order in perturbation theory, projected on the plane (r, it) and for a 
space dimension d < 4. 



*,= - + —; z = 2; 77=1; 77 = 0. (57) 

To one-loop order, the three Z factors and the parameter c a are not renormalized. Therefore, these calculations 
are insufficient to determine the fixed point value of c a as well as three of the four independent critical exponents of 
the theory. Calculations to two-loop order are necessary to obtain the full RG structure and the critical properties. 
Wilson's momentum shell RG scheme is not well suited to perform such calculations and it is technically far more 
convenient to perform these using a Callan-Symanzik RG scheme, adapted to the renormalization of critical oscillators. 
This is discussed in the next section. 



VI. CALLAN-SYMANZIK RG SCHEME 



The Callan-Symanzik RG scheme avoids the introduction of a cut-off in the momentum space, which is responsible 
for making the evaluation of multiple integrals technically difficult in Wilson's scheme. In its absence, the calculation 
of two-loop and higher order Feynman integrals is easier. Here we present the general structure of the Callan-Symanzik 
RG scheme adapted to the study of coupled oscillators. We discuss the RG flow to two-loop order in perturbation 
theory and show that the Callan-Symanzik RG scheme described here is consistent with the momentum shell procedure 
presented before. 



General formalism 



Within a Callan-Symanzik RG scheme [T5] , we start from a bare theory that follows a dynamics of the form ( 14 1, 
whose parameters are labeled with a superscript "0" . We then define the renormalized theory such that its effective 
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action is of the form (17 1. This requires to introduce a phase shift 
fields (<Pa>9a) an d the renormalized fields (<f> a ,4> a ), such that 

1/2 , 



and a frequency shift Sloq between the bare 



0°(x,t°) = {l a f,(-6uot)Zl /2 Z u <f>f,(x,t) 

4>° a (x,t°) = n a0 (-8e-5u o t)zl /2 z u 4> ( x ,t). 



(58) 



Here we have introduced Z factors for the renormalization of the fields and the time (t° = Z^t). We furthermore 
introduce dimensionless coupling constants g and g a and a scale parameter /i with u = //(47r)~ e / 2 g and u a = 

g° = Z g (Z $ Z*)-i/2g, gl = Z ga (Z^)-^g a , r" = r° + Z r r and c L 
parameters g, g a and c a on defines three beta functions. Denoting g = 

/3 = (f3,f3 a ,(3 c ) and (c^o denotes differentiation with fixed u°, u° a and c° a 
under renormalization and that we choose units such that c = 1 in the following. The renormalized correlation and 
response functions G a p = {4> a <t l i3)c and 7 a/3 = (0c* 

4>p) c are related to the physical observables C a p and Xa/3 via 



g a . Depending on /i, we relate the bare quantities to the renormalized ones by additional Z factors: 

3 ^ -1 / 2 /j_. r° = r° 4- Z„r and r° = Z Ca c a . The dependence of the renormalized 

(9, 9a, c a ), we write (3 (g, e) = ^(d M <?)o, where 
. Note again that c and D are kept constant 



The frequency ujq = Z^ 



C a/3 (x,t°) = ft Qa (-co o <)Z Z 2 G CT/3 (x,i) 

Xa/9(x,t°) = A- 1 fi Qff (0- Wo i)(Z^) 1 / 2 Z2 7CT/3 (x,<). 

5wo and the phase 9 = 6° + 58 are renormalized according to 



/u(d M (50) o 
//(9 M (5wo)o 



70(5,5a,c a ,e) 
r 7^ (g,9a,c a ,e), 



(59) 



(60) 



which defines the Wilson's functions 79 and 7 Wo . In addition, we define the Wilson's functions associated with the 
dependence of the Z factors on /i: j r — ^(9 M In Z r ) , 7^ = n(d^ In Z u )o, 7 = ^(9^ In Z^) and 7 = /x(c^ In Zt) . 

The independence of the bare theory with respect to the scale parameter fi leads to the Callan-Symanzik equations, 
which we write in terms of the renormalized theory: 



d 



din fx 



+ 



d_ 
dg 



N(N) 

»0')=i 







d 



dim 



'i(j) 



d In 5 



p(JV,iV) _ 



iV iV 
(7V + 7V-l) 7tJ + -7+-7 



p(JV,AT) 



(61) 



Here T^ N ' N ^ is the vertex function with N and iV truncated external legs corresponding to the fields 4> an d </>, 
respectively. It is a function of TV variables (wj, qi) and variables (u/j, qj) describing the frequencies and wavelengths 
associated with all external legs, and depends on the renormalized set of parameters. 



In order to calculate the Wilson's beta and gamma functions which appear in the Callan-Symanzik equations (61 1 
we decompose the bare action associated with the fields 0° and 



5° 



AO j,0 



SS 



SMop 4>a j <i>a 



(62) 



Here S-R\4> a ,4>a\ represents the action of the renormalized theory, <SMop[0a> <j>a] the action associated with the "mass 
operator" [(fra't'p], an d SS[<p a , 4> a ] combines the counter-terms. The integrals corresponding to the Feynman diagrams 
of the theory contain poles as a function of the small dimensional parameter e. The Z factors are determined such 
that the counter-terms absorb these poles and the effective action is finite. We write: 



St, 



/ { D4> a (t> a - (f> a [iu + cq 2 ] 4> a - 4> a [c Q q 2 ] \ 

ll e -(ATT)~ e - /2 g al3 (j) a (j)i3(f> 1 (f> 1 



(63) 



k 1 k 2 k 3 
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5S 



f [d (z u z^ - 1) j>J a - 4> a Ulo (z^Jz^cosse - l) 

J q,cj L \ V 

+cq 2 [z u ^j Z^Z^cos 56 - Z^^jz^Z^ Z Ca c a /c sin 56 - ij 



+c, 



iuj (zl^Z^Zprnde} 
a q 2 (z u yJz^Z^ cos 56 + Z u ^jz^Z^ c/c a sin 56 - V 



- M e (4^)- e / 2 / <j [(ZgZlcoaSB - l)g- Z ga Zlsm56g a ] 
I feifc 2 fc 3 



{[< ; 

+ [(^sa^ cos ^ - !).9a + ZgZ*sm59g] e^a^^ipj'j 



(64) 



>Mop 



| -0 Q Q yJz~^Z~^Z u (rZ r cos 50) j 



where 



9af3 



9 -9a 
9a 9 



(65) 



(66) 



Having determined the Z factors, we can calculate the beta functions. Writing 

P(g,9a,c a ,e) = -eg + (3 {4) (g,g a ,c a ) 

f3a(g,g a ,c a ,e) = -eg a + ffi (g, g a ,c a ) 
Pc(g,g a ,c a ,e) = /3( 4) (g,g a ,c a ), 

the functions (3^ do not depend on e and are given by 



(67) 



P. 



(4) 



9 B -k ^Sf 

dZ ga dZ ga 

9a g„ 9a g q 

dZl a dZ Ca 



dg ^ dg a 



dZ g 

dZ 9l 
9°- dc a 

Zc a + Ca 



(1) 



f)Z Ci 
dc a 




(68) 



where the superscript "(1)" indicates the coefficients of the poles in e _1 of the different Z factors in the matrix. Here 
Z g = Z g (Zj > Zj > )~ 1 / 2 and Z 9a = Z ga (Zj > Z^)~ 1 / 2 . The Wilson's gamma functions of the theory are independent of e. 
We have 



lr{g,g a ,c a ) = -g 



dZ; 



(1) 



dZ, 



(i) 



dg dg a 

and corresponding expressions for 7^, 7 and 7. Furthermore, 

dse^ 



(69) 



Je{g,9a,C a ) 



-9- 



dg 



9a- 



dg a 



. 1 / ds4 1] 

7„o (3, 9a, c a ) = - \-9-Qg- 



~9a 



d5u>£ ] ' 
dg a 



(70) 



The RG fixed points correspond to the values g* of g for which the beta functions are zero. The critical exponents 
of the theory are given by the fixed point values of the gamma functions. We have: 



^=(2+ 7r (r)) 1 z = 2+^(9*) 

V = 1 (<t) + 2 loJ (g*) V = 7(g*) + 2^(9*)- 



(71) 
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B. One-loop order calculations 



As in the framework of Wilson's momentum shell integration scheme, to one-loop order in perturbation theory, only 
the parameters ujq, r, g and g a are renormalized. We find j r — —2T>g and 7^ = 2Dg a , where D — AD/(Att) 2 (see 
Appendix [B] for details about the calculations). Furthermore, the beta functions associated with the renormalization 
of g and g a read 



Pa 



-eg-T> 

-eg a - V 



9 2 a - 9 2 - 2gg a c a . 2 
4 ff 



1 



\g 2 - 9 2 a + ZggaCa) - 6gg a 



(72) 



.1 + cS 

We find the same results as those discussed in Section V. Choosing c a as a parameter, the fixed points read: 

g* ~ e/5V, g* a ~ c a e/bV. (73) 



C. Two-loop order calculations and universality class 

To this order in perturbation theory, all Z factors and parameters of the theory are renormalized. Details about 
this renormalization procedure are given in the appendices: Appendix A gives graphic representations of the Feynman 
diagrams of the perturbation theory, Appendix B gives the main calculations steps and explicit expressions for the 
associated integrals, and Appendix C provides expressions for the different Wilson's functions of the theory, including 
f3 c (g,g a ,c a )- The renormalization of the propagator to second order does not affect the first-order results discussed 
above. The fixed point condition for the full problem can thus be written in the form p c (c*) = 0, where 



Pc(c a ) 



(74) 



and c* denotes the fixed point value of c a . Here, we used the fixed point values g* and g* given by Eq. (73 1. Note 
that to two-loop order in perturbation theory, this function p c (c a ) is independent of e. The full expression of the 
function p c (c a ) is given in appendix C Eq. (C2|, and its graphic representation in Fig. [2j 




a 



FIG. 2: The renormalization of the parameter c a is described by the Wilson's function p c , whose expression to two-loop order 
is given by Eq. ( |C2[ ) and which is displayed here as a function of c a . A fixed point of the theory is characterized by p c — 0. A 
single fixed point exists for c* = 0. 



A single fixed point of the theory exists with c* = and g* = 0. This fixed point is IR-stable. It is the same as 
the one of the real Ginzburg-Landau theory with 0(2) symmetry. As a consequence, the dominant critical exponents 
are the same as those known for the 0(2) dynamic model A. They are v ~ 1/2 + e/10, r\ ~ e 2 /50 and z ~ 2 + 
e 2 (61n(4/3) — 1) /50, associated with the relation z = 2 + (fj — 77) /2. Furthermore, at this fixed point, both functions 
70 and 7 Wo are equal to zero. Therefore, at the critical point, the effective phase and frequency become scale- 
invariant. Note however that since the renormalized fields differ from the physical fields by time-dependent variable 
transformations, the physical correlation and response functions differ significantly from those associated with the 
0(2) dynamic model. The asymptotic expressions of these functions are discussed is Section VII 
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The critical behaviors of the theory are characterized by the linearized RG flow in the vicinity of the fixed point. 
The linearized flow equations can be written as 

-77^- = g(^), (75) 
a In fx — 



where we have introduced the matrix 



The eigenvalues of this matrix are: 



e 

u= ( e/5 -e 2 /25X> | . (76) 
e 2 /50 



£ £ 

Wi = £ ; cj 2 = - ; w 3 = — . (77) 



In addition to u>x, which is known from the 0(2) symmetric dynamic model A, we find here two new universal 
quantities lu 2 and W3, which are specific to critical oscillators. 



D. Flow diagram of the theory to two-loop order in perturbation theory 

The RG flow of the theory within the Callan-Symanzik RG scheme is given by the variations of the three parameters 
g, g a and c a under renormalization. We display in Fig. [3]the projection of this flow on the plane (g, c a ) for two different 
space dimensions, above and below the upper critical dimension d c = 4. The first plot corresponds to d > 4 and the 



(a) g 




FIG. 3: Schematic representation of the RG flow of the theory, obtained to two-loop order in perturbation theory, and projected 
on the plane (g, c Q ). (a) Space dimension d > 4. (b) Space dimension d < 4. 

second to d < 4. 

Because the RG flow here is defined in an enlarged space, its structure differs remarkably from the one of the real 
Ginzburg-Landau theory. For d > 4, we find a line of Gaussian fixed points corresponding to g* = and g* = for 
any value of c*. These fixed points characterize the mean field universality classes of critical oscillators, which depend 
on the value of c* . Below the critical dimension d c = 4, a single fixed point exists with g* = e/52?, g* — and c* = 0. 
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Because of the existence of a whole line of Gaussian fixed points that change their stability at d = d Cl the RG flow 
has a singular structure near d = d c . As a consequence, the large scale behaviors of critical oscillators for d > d c can 
vary correspondingly to different values of the effective parameter c a . For d < d c however, the characteristic critical 
behaviors are always described by the single fixed point with c* = 0, relevant for this case. 



VII. CORRELATION AND RESPONSE FUNCTIONS AND VIOLATION OF THE 
FLUCTUATION-DISSIPATION RELATION 



In the previous section, we have discussed the RG flow and fixed point structures of the theory. We have seen 
that the renormalized fields are described in a reference frame that oscillates with rcnormalizcd frequency and phase 
factors. The fixed point theory is formally equivalent to the one of a critical point at thermodynamic equilibrium, 
namely the critical point of the dynamic Ginzburg-Landau theory with an 0(2) symmetry. However, we show now 
that the correlation and response functions of the physical fields studied here have different properties. 



A. Asymptotic behaviors of the correlation and response functions in the critical regime 



The asymptotic behaviors of the correlation and response functions of the theory near criticality can be determined 
using the RG flow and applying a matching procedure to link these functions with their expressions off criticality 
(see e.g. [38]). In the present case, such a matching procedure needs to be adapted. Indeed, the physical correlation 
and response functions are related to those defined for the renormalized fields by the time-dependent transformations 
(59 1. We therefore have to add to the usual matching procedure a scale-dependent transformation to describe the 

This transformation depends on the effective frequency and phase of 
Taking all this into account, we can write effective 



physical theory in its original reference frame, 
the oscillators, which are renormalized by the RG procedure, 
asymptotic expressions for the functions associated with physical quantities. 

The effective linear response function of the physical theory behaves as (for g£ 3> 1 and for stimulation at the 
effective frequency u)q) 



1 



1 



-n 2A r 



J8{q) 



c C ff + ii(q) 



(78) 



Here, we denote by 0Jq the effective oscillation frequency at the bifurcation, and by £ the correlation length in the 
non-oscillating phase. Furthermore, we have introduced the functions 8(q) ~ 6 e fi + a c fiq U2 + f3 c gq UJ3 and 7(g) ~ 7 e ff9^ 3 
of the wave number q — |q|, as well as non-universal effective quantities denoted by the index "eff". These functions 
are derived respectively from the renormalizations of the parameters 8 and c a in the vicinity of the fixed point. Note 
that they depend on the universal critical exponents given by Eq. ( 77 1 . Similarly to Eq. ( 78 1 , the correlation function 
behaves as 



J o ) 



1 



qZ+2-r, c 2 



cff 



i(q? 



(79) 



Related expressions can be obtained for the frequency dependence for q = in the regime (ui — uiq )£ z 3> 1. They are 
given by 



X(q = 0,w) ~ ±i 



C(q = 0,u) 



2A 



ell' 



\U — U, 



cll'l 



,eff 



DeS 

2 + z- V 



(80) 



where 



rr I & 

1 + p eS \U — Uq I 2 



(T e ff |W 
ffl 



cir 1 27, 



8(uj -u)q)~ 6» off + a cS \u - \* + Pee w - ujq 



.cff I 



(81) 



and where "±" corresponds to being positive or negative, respectively. The anomalous dependences on 

frequencies as given by Eq. ( 81 1 are due to the non-trivial evolutions of the parameters 9 and loq under renormalization. 
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B. Generalized fluctuation-dissipation relation in the critical regime 



A collective system close to a Hopf bifurcation operates far from thermodynamic equilibrium. Therefore, the 
correlation function C a /3 and the linear response function Xa/3 do not obey the fluctuation-dissipation (FD) relation 
that is characteristic of thermodynamic equilibrium. Interestingly the effective theory at the RG fixed point, expressed 
in terms of the renormalized fields <p a and <f> a , is formally equivalent to a fixed point theory at thermodynamic 
equilibrium. Therefore, exactly at the fixed point, a relation appears between the correlation and response functions 
G a p and 7 Q( g of the renormalized fields. It takes the form 



^a(3 - — 7 a /3! 



(82) 



where 7 a ^ = j' a p + *7„« has been split into its real and imaginary parts. The emergence of this relation in the critical 
regime can be discussed by the introduction of the function [53] 



2D 7i'i(q,w) 



lu Gn(q,w) 

The evolution of this quantity under renormalization is described by the following Callan-Symanzik equation: 







d In fi 



d 



d 



d 



7a 



d In uj 



7r 



d\nr 



F = 



7a> + ^(7-7) 



F. 



(83) 



(84) 



Since the fixed point theory obeys the FD relation, we have F(q, u>,r, g* , fi) — 1. Note that in mean field theory we 
find F mf = (lu 2 + R 2 - (c a q 2 ) 2 ) / [uj 2 + R 2 + (c a q 2 ) 2 ), which differs from F = 1 if c a ^ 0. 

The fact that the renormalized theory at the fixed point obeys the FD relation F — 1 , implies that the correlation 
and response functions C a p and Xap &t that point are not independent. Since they are related to G a p and j a g by 



Eq. (591, we find 



cos 6> cff xi'i + sin 9 c sXi2 = 



cos0 e ffXi2 _ smfcffXn 



1 



1 

2A c ffD c ff 



(ufCu + iloC 12 ) 



(85) 



Here again \afi — x' a p + i-x'Lp nas been separated in its real and imaginary parts. At the bifurcation, this relation is 
asymptotically satisfied in the long time and wave-length limits. It is a consequence of symmetry properties of the 
fixed point theory which impose constraints on the correlation and response functions at criticality. Indeed, the FD 
relation is connected with time-reversal invariance, which emerges for the fields 4> a and <p a at criticality while it is 
not obeyed for the physical fields ip a and ip a . 



C. Breaking of the fluctuation-dissipation relation 



The relation between the physical correlation and response functions at criticality given by Eq. (85 1 is not an 
FD relation. In order to characterize the violation of the FD relation between C a p and Xap, we define an effective 
temperature T g, which depends on frequency and wave vector [39] : 

Jrff(q,^) = u Cn(q,m) 

Here, kg denotes the Boltzmann constant and T is the temperature. Using the previous asymptotic expressions for the 
two-point correlation and response functions, we find universal behaviors of this effective temperature at criticality: 

T cS ( qi LO = iof)/T ~ q- z 

T cS (q = 0,u;)/T ~ \uj-ujf\~ a . (87) 

For the particular case c a — and u a — 0, a — 1, while otherwise a ~ 1 — e/5 to first order in e|54]. This singular 
behavior of the effective temperature implies a violent breaking of the FD relation. This is consistent with the fact 
that spontaneously oscillating systems operate far from thermodynamic equilibrium. 
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VIII. SUMMARY AND CONCLUSION 



We have studied the critical behaviors of a large number of locally coupled oscillators when approaching a homo- 
geneous synchronization transition from the disordered phase in a d = 4 — e dimensional space. On large length 
and time scales, the critical behaviors can be described by a statistical field theory that is given by the complex 
Ginzburg-Landau equation with additional noise and forcing terms. At the critical point of a homogeneous oscillatory 
instability, time-translational invariance is spontaneously broken in the system. The field variable Z in the complex 
Ginzburg-Landau field theory is constructed in such a way that time-translations correspond to global phase changes 
of this complex variable Z. Within this framework, breaking of time-translational invariance becomes formally similar 
to the traditional spontaneous symmetry breaking known for other second-order phase transitions. 

We have established the structure of the associated dynamic RG within Wilson and Callan-Symanzik schemes, 
and performed the calculations to two-loop order in perturbation theory. We have shown that the critical point is 
formally related to the equilibrium phase transition in the real Ginzburg-Landau 0(2) dynamic model A. However, 
the RG flow of critical oscillators is defined in a larger parameter space of non-equilibrium field theories and leads 
to a renormalization of oscillation frequency and phase. The FD relation is broken in the system, which can be 
characterized by an effective frequency-dependent temperature, diverging at the effective oscillation frequency with 
an anomalous power-law. 

The formal analogy with an 0(2) symmetric dynamic field theory, valid at the critical point, leads to several 
interesting results. For d > 2, the collective dynamics of coupled oscillators exhibits a second order non-equilibrium 
phase transition. This phase transition is a generalization of Hopf bifurcations, which are conventionally defined in 
the context of nonlinear dynamics, to non-equilibrium statistical physics. On the oscillating side of the bifurcation 
and in the thermodynamic limit, the system exhibits long-range phase order and coherent oscillations. In mean field 
theory, the universal properties of this oscillating instability are captured by the normal form known from nonlinear 
dynamics. Below the upper critical dimension however, fluctuations become relevant and anomalous scaling laws and 
critical exponents appear. 

The case of critical coupled oscillators studied here provides a further example for the emergence of an equilibrium 
universality class in a non-equilibrium dynamic field theory. In non-equilibrium systems with non-conserved order 
parameter, detailed balance is often effectively restored at criticality [T7]. This is the case e.g. for the model A 
dynamics of the real Ginzburg-Landau theory with Z<i symmetry (401 141) , even when the symmetry is broken by the 
non-equilibrium perturbations |42j . and for some of its generalizations to the 0{n) symmetry |43j . In the present 
case, the detailed balance condition is not restored for the physical variables, but appears only in the oscillating 
reference frame associated with the effective frequency and phase of the oscillations at the transition. This emergence 
of detailed-balance symmetry at criticality imposes a generic relation between the correlation and response functions 
of coupled oscillators as given by Eq. (85 1. 

The structure of the RG flow studied here is singular at the upper critical dimension d c = 4. Indeed, as depicted 
in Fig. [3j the line of Gaussian fixed points, which is stable above d = 4, becomes unstable for d < 4 where only one 
isolated stable fixed point remains. Our results obtained in an epsilon expansion are valid close to the upper critical 
dimension for d = 4 — e. We can speculate how our results are modified in lower dimensions d. In analogy with 
the equilibrium 0(2) dynamic model, we expect the phase order of the oscillations to vanish for d < 2, and to be 
quasi-long range exactly at the lower critical dimension d — 2. In the last case, spectral peaks on the oscillating side 
of the Hopf bifurcation are expected to exhibit power-law tails with non-universal exponents. If the formal analogy 
with the equilibrium critical point found here in d = 4 — e persists in d = 2, we would expect to see features of the 
Kosterlitz-Thouless universality class |H] in systems of coupled oscillators in this dimension. 

The different values of the space dimension d of coupled oscillators can be related to different realizations of coupled 
nonlinear oscillators in various physical and biological systems. The mean field limit d > 4 is found in systems where 
oscillators are coupled by long-range interactions. Examples for such a situation are sarcomere oscillations in muscles. 
There, large numbers of myosin motor proteins generate oscillations when interacting with actin filaments, which 
represent tracks along which the motor proteins move. Oscillations occur if the motor collection acts against elastic 
elements and in the presence of a chemical fuel that supplies the necessary energy. For stiff filaments, this situation 
is well described by globally coupled motor proteins for which mean field theory applies at the Hopf bifurcation [IS] . 

Systems of coupled oscillators in three dimensions could be realized in oscillatory chemical processes in bulk solution. 
This is the case e.g. of the Belousov-Zhabotinsky reaction which can be studied in the framework of the complex 
Ginzburg-Landau equation (see e.g. [46J. On mesoscopic scales, the system can be viewed as a collection of interacting 
volume elements, each representing an individual chemical oscillator. There, oscillations are subject to fluctuations 
due to the finite number of reacting molecules present in each volume clement . 

Coupled oscillators in two dimensions can be realized by oscillators arranged on a surface. Such a situation may 
occur in the clcctrosensory organ of some fish species where many electrically oscillating cells constitute the sensory 
epithelium [47, 48J. Critical oscillators coupled in two dimensions can also in principle be realized in artificial systems. 
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Nanotechnology aims to build functional units on the sub-micrometer scale. Large arrays of nano-oscillators on 
patterned substrates coupled to their neighbors by elastic or viscous effects would provide a 2-dimensional realization 
of our field theory. Finally, the case d — corresponds to a single noisy oscillator. Here, fluctuations destroy the 
Hopf bifurcation and only its signatures can be observed. In the context of biological systems, an example is the 
spontaneous oscillations of the mechano-sensory organelle of auditory hair cells 39, 49 . Here, the critical divergence 
of the linear response function is ideally suited for signal detection. 

In order to observe the critical exponents discussed here, homogeneous chemical oscillations in a bulk system with 
d — 3 (which corresponds here to e = 1) would be a good candidate. However, the critical exponents attached to the 
RG fixed point are only observable when the system is observed sufficiently close to the critical point. The range and 
experimental accessibility of this critical regime can be estimated by a Ginzburg criterion, see Appendix [D] Assuming 
that in a chemical system, the Hopf bifurcation occurs if a molecular concentration p exceeds a critical value p c , the 
critical regime corresponds to 

\P~ Pc\ W 



p c p 2 T 2 C 3 ' 

where u>o is the oscillation frequency and t c denotes a chemical reaction time. The coefficient c here is the bare 
coefficient describing the coupling of oscillators in the complex Ginzburg-Landau equation. We can rewrite this 
expression as 

J ^<7^(— Y> (89) 



(Pl 3 c) 2 



C 



where I 2 — D m r c is a reaction length and D m denotes a microscopic diffusion coefficient. Since ojqt c <C 1 (oscillations 
are slow compared to fast reaction times) and l\p 3> 1 (the volume per molecule is small compared to the reaction 
volume), accessibility to the critical regime requires that D m ^> c. This condition is satisfies if c becomes small. 
This happens in particular if a Turing instability is approached. At the point where such an instability occurs, the 
coefficient c changes sign and stationary spatial patterns appear. Our analysis suggests that before such a point is 
reached, the critical regime of the Hopf bifurcation becomes accessible. Therefore, the scaling behaviors and critical 
exponents discussed here could be experimentally observable in oscillating chemical systems. 
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APPENDIX A: FEYNMAN DIAGRAMS OF THE PERTURBATION THEORY 

Here, we present the graphic representation of the terms of the expansion scries that we used for the calculations to 
one and two-loop orders in perturbation theory. The expansion series of the action of the complex Ginzburg-Landau 



theory given by Eq. ( 17 1 can be represented by Feynman diagrams as usual. The free propagator of the theory, 



calculated from the Gaussian part of the action of the decomposition (62 1 in a Callan-Symanzik RG scheme, is 
given by 



2D ( oj 2 + (c Q q 2 ) 2 + (cq 2 ) 2 2icjc a q 2 



A| 2 V -2iujc a q 2 u 2 + (c Q q 2 ) 2 + (cq 2 ) 2 



n r \ 1 ( —iuu + cq 2 c a q 2 \ , 
^> W ) = a{ -cW -J"cq 2 j' < A1 > 

where A = (— iu> + cq 2 ) 2 + (c a q 2 ) 2 . The interaction vertex reads: 

- U a p6 j5 = -^(4TT)- e/2 g af3 S l5 , (A2) 

where g a p is given by Eq. (66 1. Graphic representations of these elements are displayed in Fig. [i] Due to the 
presence of the non-diagonal clement u a in the matrix U a p, the interaction vertex contains three non-equivalent types 
of external "legs". The symbol used for the interaction vertex in Fig. [3] indicates this fact. The expressions of the 



counter-terms of the theory are directly visible on the decomposed expressions ( 64 ) and ( 65 1 of the action of the 
theory. 
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FIG. 4: Graphic representation of the propagators G a p and 7a/3, and of the vertex U a p5 1 s. 



Fig. 5 displays the diagrams that contribute to the renormalization of the "mass operator" \<t> a (f)p\ and the vertex 
to one-loop order, and Fig. [6]shows the diagrams that contribute to the renormalization of the propagator to two-loop 
order in perturbation theory. 

(a) 




(b) 





FIG. 5: Feynman diagrams contributing to the renormalization of the theory to one-loop order, (a) Renormalization of the 
"mass operator" [<f> a cj>p\. (b) Renormalization of the interaction vertex. 



(a) 
(b) 

WW(A^ — vw^^ — wwCi^ 

VWfi^ Wrt^^ 

FIG. 6: Feynman diagrams contributing to the renormalization of the propagator to two-loop order in perturbation theory. 

(a) Contribution to the renormalization of the noise amplitude parameter D(Z UJ Zj > — 1) in the counter-terms SS in Eq. (64 1. 

(b) Contribution to the renormalization of the other counter-terms to the propagator. 



APPENDIX B: EXPLICIT EXPRESSIONS OF THE CALCULATED PERTURBATION SERIES 

We present here explicit expressions of the Feynman integrals associated with the diagrams displayed previously, 
to one and two-loop order in perturbation theory. Using the notation introduced in Subsection |VI A[ we denote 

ky ^{I'^'ip -} th° contribution of the n-loop order to the vertex function with N and N truncated external legs 

corresponding respectively to the fields <f> ai and 0^., and with L insertions of the mass operator [<f> a (j)p]. 
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1. One-loop order 



To one- loop order in perturbation theory, the propagator of the theory is not renormalized. The renormalization 
of the "mass operator" [4> a (j)p\ is given by 



r (i i,i)(i) = r ^£( 47r ) e/2 Jd [ g ^s aa + gf3a 8 ai + gp a 8 1CT ] 5 1<T , 



(Bl) 



where g a p is given by Eq. (66 1. The expression of the integral J d is given below. Because of phase invariance 
symmetry, described for p^'W') by Eq. (jfjj, only two terms need to be calculated. 

The renormalization of the vertex function to one-loop order in perturbation theory at a symmetry point of the 
configuration of external momenta, respects the following symmetry: 

rSjS 3 (fci,fe,fe, k 2 ) = [M ai01 6^ 3 + Af Ql/32 <Wa + M ai03 8 M ^\ x (27r) d+1 ^+ 1 (fc 1 + 3fc 2 ), (B2) 

where M a p has the phase invariance symmetry (16JI. Therefore, the renormalized interaction vertex of the theory 
remains of the same structure as the original one, and only two independent terms need to be calculated, e.g. 



r 



(1,3,0)(1) 

l.iii 



and r 



(1,3,0)(1) 
2.111 



For vanishing external frequencies, we get 



1 



1 



9 ~h + iJd)+gA-7ih)+ 2gg a ~i d + u A \+ c.c 



I d + 2iJ d ) + 2gg a [ -I d + J d ) + c.c 



ff -Tih 



where "c.c." denotes the complex conjugated value. 

In the previous expressions, the integrals I d and J d are given by 



I d (c,c a ,q) = j 
Jp 

J d {c,c ai q) = 

J p 



D 



cp 2 (c + ic a ) [p 2 + (p - q) 2 ] 
D 

cp 2 [(c - ic a )p 2 + (c + ic a )(p - q) 2 ] 
Within a dimensional regularization scheme, they read: 



„d-4 



Id(c,c a ,q) = 



Jd(c, c a ,q) 



r 



4-d 



f ^_ 

c(c + ic a ) (4tt)<V 2 

D 1 1 

c(c + ic a ) (47r) 2 (47r)- £ / 2 ~e 

D q d - 4 (d\ {A-d 
c (47r) d / 2 ^2 



x 2 



o (l + z) rf ~ 2 



fi.T 



(l + 0(e)) 



x 



1 dx 
(1 - x + 2cx) 2 



1 - x + x(c-ic a ) 
X{c+Ka) (l-x + 2cx)* 



1 



1 



c 2 (47r) 2 (47r)- e / 2 



where q = |q| 



(B3) 



(B4) 



(B5) 



2. Explicit expressions of the Feynman integrals to two-loop order 

Since we are looking for the first non-trivial corrections to the critical behaviors in perturbation theory, we only 
need here to renormalize the propagator, which to one-loop order was not renormalized. Following the same notations 
as previously, and for vanishing external frequencies, we have 

K2) (k) = M 2e (47r)- V + g 2 a ) [21 A + 6I B ] S a0 

r (M,0)(2) (k) = 2(U 2 £(47r) - £ ^2 (4Ja + 7Jc + j e) + g 2^j A _ 7Jc _ j e) + 2 gg a (7J D - J F )} 

rgi 1,0)(1) (k, 0) = 2 M 2£ (4tt)- £ [g 2 (-AJ B - 7J D + J F ) + g 2 a {-AJ B + 7J D - J F ) + 2gg a (7J c + J E )] . (B6) 
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Here 



I A = 




(c Q ) + id.e'e" 


-Ca)) 


e' = 


+1 


e" 


= -1 


Ib = 




(c Q ) + 7<j, e ' e " 


-Ca)) 


e' = 


-1 


e" 


= +1 


J A = 




>{c a ) + Jd,t't" 


(-Ca)) 


e = 


+1 


e 


= -1 


Jb = 


-^(Jd,e'e' 


'(c a ) - Jd,e'e" 


(~Ca)) 


e = 


+1 


e 


= -1 


Jc = 


^{Jd,e'e- 


'(c a ) + Jd,e'e" 


(-Ca)) 


e = 


-1 


c' 


= +1 


Jd = 




'(c a ) - Jd,e'e" 


(-Ca)) 


e = 


-1 


e 


= +1 


Je = 




'(c a ) + Jd,e'e" 


(-Ca)) 


e = 


-1 


e 


= -1 


Jf = 


^(Jd,e'e' 


>(c a ) - Jd,t't" 


(~Ca)) 


e = 


-1 


e 


= -1, 



(B7) 



and 



I d ,e'e"(c, C a ,k) = / 

"'p.q 

Jd.e'e"(c, C a ,k,Wfc) = / 



D 3 



cq 2 c(p - q) 2 c(k - p) 2 [(c + ic a )q 2 + (c + ie'c a )(p - q) 2 + (c + ie"c a )(k - p) 2 ] 

D 2 

cq 2 c(p - q) 2 [(c + ic a )q 2 + (c + ie'c a )(p - q) 2 + (c + ie"c a )(k - p) 2 + iu k ] ' 



(B8) 



The expressions of these integrals as a function of the space dimension d are too large to be displayed here. We 
therefore just report the expressions of their divergent parts as e goes to zero: 



D 



1 



1 



c 4 (47r) 4 (47r)- e e 



I d ,e'e"(c, C ,k) 

dJ d ,e'e"(c, c a ,k,o; fe ) 

<9ck 2 ~~ c 4 (47r) 4 (47r)- £ e 

dJ d ,eU"(c, c a ,k,uJk) D 2 1 1 



ff eV ,(c Q )(l + 0(e)) 



Z? 2 



1 



1 



ifl,eV'(Ca)(l + 0(e)) 



c 3 (47r) 4 (47r)- £ e 



(B9) 



where c a = c a /c. The integrals i? e ' £ "(c a ), ifi ie 'e"(c ) and i^2,e'e"(c Q ) are given by expressions that are similar to the 
ones found to the one-loop order (Eq. (B5l), which for d = 4 reduce to integrals of rational fractions with complex 
parameters. The result of these integrations read: 



Kl t + l + l(c a ) 


= h(Ca) 


Ki-x-i(c a ) 


= h(Ca) 


Ki (Co) 


= h(Ca) 


Ki t+ i-i(c a ) 


= h(c a ) 


K 2 , + l + l(c a ) 


= h(Ca) 


K 2 ,-l-l(Ca) 


= h(Ca) 


K 2 -l + l(Ca) 


= h(Ca) 


K 2 .+i-i(c a ) 


= h(c a ) 


H+l-l(5 ) 


= #-1+1 



1 


6(1 


+ ic a ) 


2- 


f ICa 


4(3 


+ ICa) 


h(~ 


-C a ) 


1 - 


ICa 


6 - 


2ic a 




1 


(H 


" ICa) 2 




1 


(1- 


" iCa) 2 


h(- 


-C a ) 




1 


(1- 


- iCa) 2 



3 + ic a 



(3 - ic a )(l + ic ) 



(BIO) 
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where Ig(c a ) can be represented by 
h{c a ) 



(1 + 3) 



4 In [2] — i arctan 
In [-i (1 + 2ic a +cl)] + In 
1 



2Cg 

3 + 52 



~la[9 + 10c£ + c£] 



. 1 + 2ic a + C 2 a 



lC a 



In 



3 + 2ic a + c* 

y/{l + ic a ) 2 (l-ic a ) 
l + ic a + + ic a )'- 



In 



, 2 + ic a + cl- + ic a f 
'2 + ic a + c 2 a + ^(1 + iCa) 2 



(Bll) 



Because of the ambiguity in the definition of the complex logarithm, the integrals Iq{co), It{co), I%(c a ) and Ig(c a ) are 
not uniquely defined by these expressions [55 . To specify entirely these functions, one needs to use their values for 
c a = 0, which read: 



h{c a = 0) = I 7 {c a = 0) = h{c a = 0) = In 
h(c a =0) = 3 In 



(B12) 



The entire functions are then defined as the unique analytic prolongations of the expressions (BIO I and (Bll ), defined 
in the vicinity of c a — together with the specification (B12l. 



APPENDIX C: EXPLICIT EXPRESSIONS OF THE RENORMALIZATION GROUP EQUATIONS 

With units such that c = 1, the function /3 c (<jf, g a , c a ) is given by: 

(3 c (g,g a ,c a ) = V 2 [(g 2 + g 2 a )lm(-h) + (g 2 - g 2 a )lm(2I 2 ) + 2gg a Re(-2I 2 ) 
+V 2 c a [{g 2 + g 2 a )Re(h) + (g 2 - g 2 a )Re(2I 2 ) + 2gg a Im(2I 2 ) 
+V 2 c a (c a + c" 1 ) \(g 2 + g 2 a )lm{h) + {g 2 - g 2 a )lm(-2I e ) + 2 5ffQ Re(2/ 6 )l , 



(CI) 



where T> = 4Z?/(47r) 2 [55]. The function p c (c a ), which is defined in Eq. (74 1 and determines the fixed points of the 
theory to too-loop order, is given by: 



Pc(c a ) 



50 



(l + cgW4(l-cg)Im 



In 



\3+tc a J 



{1-iCaf 



2(1 + cl)lm 



In 



( i \ 

(1 - iCa) 2 









— 8c a Re 




{l-ICa) 2 





(C2) 



where the three complex logarithms are defined by prolongation of their real values for c a = 0. This function is 
displayed in Fig. 

Finally, the other Wilson's functions of the theory are given by the following expressions: 



le (g,g a ,c a ) =V 2 (g 2 + g 2 a )Im(I 8 ) + (g 2 -g 2 a )Im(-2I 6 ) + 2gg a Rc(2I 6 ) , 



(C3) 
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and: 



l(g, 9a, c a ) = V 2 (g 2 + g 2 a )Re(3h - h - h) + (g 2 - gl)R^h - 2/ 6 ) + 2 55Q Im(6/ 2 - 2/ 6 ) 



+32? 2 c a 



(g 2 + g 2 a )lm(I s ) + (g 2 - s =)Im(-2J 6 ) + 2gg a Re(2I 6 ) 



M9, ga,c a ) = V 2 y{g 2 + g 2 a )Rc{-I 4 + h) + {g 2 - g 2 a )Re(-2I 2 + 2I 6 ) + 2. 9 .gJm(-2/ 2 + 2/ 6 ) 
-V 2 c a \{g 2 + gl)lm(I s ) + {g 2 - 5 2 )Im(-2/ 6 ) + 2gg a Re(2I 6 ) 



7(3, g a ,c a ) = V 2 [(g 2 + g 2 a )Rc{h ~h + h) + (g 2 ~ 9l)M^h - 2J 6 ) + 2gg a lm(2I 2 - 27 6 ) 
+V 2 c a \{g 2 + &)lm(h) + {g 2 - g 2 a )lm(-2I 6 ) + 2gg a Re(2I 6 ) 



(C4) 



APPENDIX D: GINZBURG CRITERION 



The critical behaviors described by the RG fixed point are valid in the proximity of the critical point. In practice, one 
can estimate how close to the critical point an experiment needs to be performed in order for the critical behaviors to 
become observable. Further away from the critical point, when mean field theory is still an appropriate approximation, 
nontrivial critical exponents are unobservable. The Ginzburg criterion estimates the breakdown of mean field theory 
at the point where the variance of order parameter fluctuations exceeds its average |50j . For simplicity, we discuss 
here the real equation ^ with u a = and c a = 0. Based on this criterion, mean field theory is valid if 



-(d-2) 



< 



(Dl) 



with £ 2 = c/\r\. We now relate this expression to a bulk chemical system where oscillations occur if the concentration 
p of some species exceeds a critical value p c : r ~ a(p c — p) where a is a proportionality coefficient. In this case, the real 
part of Z is related to the concentration fluctuations Sp/ p. As a consequence, r and u have dimensions of an inverse 
time, while the noise strength D is a volume per unit of time. The coefficient c has units of a diffusion coefficient. We 
can estimate the fluctuations of the number N of molecules in a reference volume v after a time t using dtSp ~ prj as: 



(5N 2 



d d xd d x'(Sp(x)Sp(x')) ~ p 2 Dvr. 



(D2) 



We estimate (SN 2 )/N 2 ~ l/N ~ within a reaction volume v = 1% and a reaction time r c , with I 2 — D m r c and 

where D m is a microscopic diffusion coefficient. From this, it follows that D ~ l/(pr c ). For d = 3, Eq. (Dl I implies 
that mean field theory is valid if 



\P-Pc\ u 2 D 



r c 



3 ' 



(D3) 



where rp = ap c . Estimating ro ^ u ^ ujq, we find Eq. (88 1 and (89 1 
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